Macroscopic Kinetic Effect of Cell-to-Cell Variation in Biochemical Reactions 
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Genetically identical cells under the same environmental conditions can show strong variations 
in protein copy numbers due to inherently stochastic events in individual cells. We here develop a 
theoretical framework to address how variations in enzyme abundance affect the collective kinetics 
of metabolic reactions observed within a population of cells. Kinetic parameters measured at the 
cell population level are shown to be systematically deviated from those of single cells, even within 
populations of homogeneous parameters. Because of these considerations, Michaelis-Menten kinetics 
can even be inappropriate to apply at the population level. Our findings elucidate a novel origin of 
discrepancy between in vivo and in vitro kinetics, and offer potential utility for analysis of single-cell 
metabolomic data. 

PACS numbers: 87.18.-h, 82.39.-k, 87.16.-b 



Noisy or stochastic molecular events are plentiful in 
the life of a cell. In the past several years, extensive 
experimental and theoretical efforts have been devoted 
to analyzing stochastic processes in gene and protein ex- 
pression, at both the transcriptional and translational 
levels [l[ . One consequence of stochastic gene expression 
is that the number of molecules of a given protein can 
vary substantially from cell to cell, even within geneti- 
cally identical populations @, Q. Such stochastic gene 
expression has received considerable attention in relation 
to cellular regulation, phenotypic diversity, and disease 

gt 



Kinetic modeling of metabolism aims at achieving a 
quantitative description of biochemical reactions to gen- 
erate mass and energy required for cell survival. Since 
there are rarely available techniques to detect metabo- 
lites inside a single cell most kinetic models for in- 
tracellular metabolism have been built on experimental 
data obtained from cell populations. These models, how- 
ever, would be valid when behavior of individual cells is 
very similar to the average behavior of the population. 
Significant variation of enzyme abundance between cells 
might challenge this traditional approach, and here we 
investigate how single-cell variation affects the kinetics 
of metabolic reactions appearing at the population level. 

The flux of any given reaction (y) can be expressed 
as a function of the molecular concentrations and the 
kinetic constants associated with the reaction: v = 
f(E, {5 7 }, {Ki}), where E stands for the concentration 
of enzyme, Sj for the concentration of the 7th metabo- 
lite participating in or allosterically regulating the re- 
action, and Ki for the ith kinetic constant in the re- 
action. For example, if a single-substrate reaction fol- 
lows the Michaelis-Menten (MM) kinetics, then v = 
f(E, S, {K , K M }) = K ES/(K M + S). Usually, exper- 
iments measure concentrations or fluxes averaged over a 



cell population, as appears below. 

(v z ) z = (f(E z ,{S lz },{K t })) z 

= f({E z )„{(S^),},{K i })x(l + t), (1) 

where we index each cell with z such that v z is the re- 
action flux inside cell z, and (• • • ) z averages a value over 
a cell population. The last equality gives the definition 
of £ to be equal to zero if molecules are uniformly dis- 
tributed over the population. For clarity of analysis, here 
we only consider cellular variability of molecular concen- 
trations but not of kinetic parameters and cell vol- 
umes; in Eq. (JXJ) , there is no z-dependency of Ki and 
no weighting of cell volumes for the averages. Although 
this analysis can be extended straightforwardly for the 
factors excluded here, our work demonstrates that ex- 
perimentally well-established heterogeneity in single-cell 
enzyme abundance itself gives rise to inevitable and 
definite effects on the kinetic properties of metabolic re- 
actions. In the following, we consider the steady states of 
reactions under a given enzyme distribution, as the typi- 
cal time scale of reaction rates is much shorter (~ 10~ 3 s) 
than that of enzyme concentration changes (> 10 2 s) [8|. 
Applying the series expansion, we can write £ as 
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where the higher-order terms have been omitted, fo = 

/ = /(£,{S 7 },{i^}), AE Z = 
E z — (E z ) z , AS JZ = S yz — (S JZ ) Z , and |o is for derivative 
around E = (E z ) z , S 7 = (S JZ ) Z . Note that (i) if / is 
a simple linear function of molecular concentrations, £ 
vanishes, and (ii) variability of the concentrations across 
cells explicitly contributes to £. In other words, inherent 
nonlincarity in reaction kinetics and significant cellular 



2 



variation in molecular levels manifest the difference be- 
tween {v z ) x and f((E z ) z ,{(S lz ) z }, {Kj}). Hence, kinetic 
constants calculated from a population do not necessar- 
ily coincide with {K{} at a single-cell level, as will be 
addressed in detail. 
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Suppose that {v z ) z , 
perimentally measured 

stants {Ki} fitting f((E z ) z , {<5 7f ),}, {Ki}) 

for example, by minimizing 

f((E z ) z ,{(S, z ) z } 7 {K t })-(v z ) z )/(v : 
the typical approach. If more than one experimental 
dataset are available, {K%\ may instead be obtained 
by minimizing J2 r e r> where subscript r represents 
the rth experiment. Combined with Eq. {jTJ) , this 
procedure results in the following formula for {Ki} at 
the lowest-order approximation: 
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where AK and u are the vectors whose elements are 



AKi = Ki - Ki 
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respectively, and A 1 is the inverse matrix of A with 
elements 
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It should be noticed that Eq. J3J gives a degree of dis- 
crepancy between kinetic constants at a population level 
{Ki} and those at a single-cell level {Ki}, as a function 
of a degree of cellular heterogeneity in Eq. ((2J . 

Taking into account a single-substrate reaction gov- 
erned by the MM kinetics [/ =JC ES/(K M + S)}, it is 
straightforward to get K and Km from Eq. (|3j) if both 
come out of the same population data: 
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where cov r (x r , y r 



{x r y r 



(x r ) r (y r ) r , var r (x r ) 



{x^) r — (x r )r- From Eq. one can further prove the 
following relation: 
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which shows that the degree of deviations in kinetic pa- 
rameters is essentially determined by (£r) r , as the right 



side ranges from |(£ r ) r |/\/2~ to |(£ r ) r |. On the other hand, 
if Km is known at the single-cell level, one might calcu- 
late only K a from experimental data, but not Km- In 
this case, the following simplified relation from Eq. §3§ 
holds: 
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Since £ r s consistently include the negative-sign terms (as 
shown below) , they would not be simply canceled out by 
each other under averaging, thereby allowing for signifi- 
cant nonzero (£ r ) r in Eqs. (|5|) and ©• For simplicity of 
analysis, we will concentrate upon cases of Eq. © out of 
single experiments. Note that S z = Kmv z /(K E z — u z ); 
thus, AS Z in Eq. © can be substituted for by AE Z and 
Av z = v z - (v z ) z to give, 
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As long as there exists variation in enzyme concentration, 
the first term on the right side of Eq. ([7J has a nonzero 
magnitude, always greater than {{AE Z ) 2 ) Z /{E Z ) 2 of 
which experimental values are recently available for the 
yeast Saccharomyces cerevisiae Q. Let £a be this first 
term, and we can estimate £ to be £a if ({Av z ) 2 ) z is 
sufficiently small. Without the lowest-order approxima- 
tion from Eq. ([2]), we can also get £ from the exact 
formula by setting ({Av z ) 2 ) z = 0, provided that the en- 
zyme concentration follows the log-normal (= £r,) or nor- 
mal (= £at) distribution 0, Sj- Such £l(n) satisfies the 
following equality: 



(S z 



K M 



P L(N){x) 



l+<5 



1 



dx . 



(8) 



where P L {N){x) = ((v z ) z /K ) x P L(N) (E), x = 
(K /{v z ) z ) x E, (v z ) z = K (E Z ) Z (S Z ) Z /(K M + (S z ) z ) x 
(1 + an d Pl(n)(E) is the probability distribu- 

tion of enzyme concentration E almost following the 
log-normal (normal) distribution. Because Pum(E) 
is approximated as the log-normal (normal) distribu- 
tion, Plin)( x ) can be approximated as the same, with 
(x ) = (l + K M /(S z ) z )/(l+£ L{N) ) and ^/{(AxJ)/{x) = 
y/({AE z ) 2 ) z /(E z ) z where Ax = x - (x). S > is cho- 
sen small enough to satisfy Jq^^s PL(N)( x )dx <C 1. Al- 
though £l(N) c& n be accurate for the particular forms of 
enzyme distribution following the log-normal (normal) 
distribution and ^a is just an estimation, £a is poten- 
tially useful as can be applied without knowing a specific 
functional form of enzyme distribution. 

We now turn our attention to the empirical values of 
£a, and Table Q] shows results for several irre- 
versible reactions, which follow MM kinetics and have the 
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necessary experimental data for calculating £4, £5, and 
£at for S. cerevisiae. It is then observed that £4, and 
£iv easily reach ~ 0(— lO^ 1 ), thereby lowering K Q by sev- 
eral ten percent of K according to Eq. ©. Furthermore, 
a striking range of the metabolite concentration changes 
(~ 100-fold increased or decreased) against different en- 
vironmental conditions Q can overweight (S z ) z /Km in 
Eq. for certain conditions, which could drastically 
increase the magnitude of £4 as well as of £. 

When can £4 be used to safely infer £? In a directed 
linear pathway as shown in Fig. QJa), each reaction takes 
as a substrate the product of the preceding reaction, and 
frequently involves an additional cosubstrate (such as a 
water molecule) that is abundant in the cell and whose 
variations can be neglected. For any given reaction in 
such a pathway, v z at a steady state is entirely deter- 
mined by the influx from the upstream region indepen- 
dently of the given reaction itself (Sv z /8E Z =0), if there 
is absent any regulatory connection between the reaction 
and the upstream region. Thus the last term on the right 
side of Eq. ((7j) vanishes because of decoupled v z and E z , 
and a negative sign of the remaining second term even 
weights the effect of £a on £ because £a also has a neg- 
ative sign. Therefore, the lower bound of the magnitude 
of £ in Eq. can be predicted by £4- 

If we consider more elaborate pathways than simple 
linear pathways, the last term on the right side of Eq. ([7]) 
may not simply vanish, and could play a role in determin- 
ing £. Specifically, if E z and v z are positively correlated, 
the last term will weaken the effect of £4 on £, and if 
they are negatively correlated, will weight the effect in 
an opposite way. Here we focus on the former case, as 
examplified in Figs. QJb) andQJc). Figure QJb) depicts 
a negative-feedback case where the substrate of a given 
reaction inhibits the first reaction in the pathway or the 
transport of its precursor. In such a way, flux is reduced 
when the substrate is accumulated, as characterized by 
v z — c z /{\ + (S z /Kj) h } where c z is the maximal influx, 
Kj is the dissociation constant of the inhibiting interac- 
tion, and h is a Hill coefficient. Substituting this into Eq. 
(O gives rise to, 



(S z 



K M 



■(1-a) 



((A£ z ) 
(Ez) 



(9) 



where a = fl + {1 + (K > / '{S z ) z ) h }h,- 1 (1 + 
(S Z ) Z /K M )~ 1 }- 1 0. Compared with Eq. Q, 
the contribution of the term with ((AE v 



l )*/(E z 



is 

relatively small, and even approaches zero at a strong 
inhibition limit {(S z ) z Kj, h ^> 1). This effect origi- 
nates from the presence of a positive correlation between 
E z and v z , as anticipated above. A similar effect can 
also be found from the case of Fig. [ljc). In branching 
pathways, such as depicted in Fig. (He), different 
enzymes can bind to the common substrate, and each 



catalyzes a first reaction in a different pathway. More 
specifically, a metabolite S z with influx c z is converted to 
a product with rate v\ = K^E^S Z /(K M + S z ) or to the 
other product with rate v 1 ^ 
Equation ([7J then leads to, 
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As expected, the contribution of the term with 
((AE I z ) 2 ) z /(E I z )l is smaller than that with 
((AE Z ) 2 ) Z / (E Z ) 2 Z in Eq. (J7J, and even approaches 
zero when (^) 2 -C (v I z I ) z - Interestingly, the second 
term on the right side of Eq. (|10|) explicitly comes from a 
correlation between E z and E* z . It should be noted that 
such a correlation between different enzymes does not 
only affect £s for specific branching pathways, but also 
does so for many general situations, because the second 
and third terms on the right side of Eq. ([7]) contain a 
variation of flux, and flux may change in response to 
the activity of enzymes or proteins other than only that 
of an enzyme of interest. Hence, a correlation between 
different enzymes as well as a variation of each enzyme 
can play an important role to distinguish a gap between 
single-cell-level and population-level kinetic constants. 



(a) 



(b) 



(c) 

< 



FIG. 1: Various types of metabolic pathways. Squares for 
metabolites and circles for reactions. Fluxes run into arrowed 
directions, (a) Linear, (b) linear with feedback inhibition de- 
picted by a dashed line, (c) branching pathways. 



It is also interesting to note that in a limited range our 
results may be experimentally accessible without mea- 
suring metabolites inside a single cell. Briefly, let us con- 
sider an enzyme taking as a substrate a molecule which 
has just passed from the extracellular medium along a 
membrane transporter like a permease. If v z cc T z where 
T z denotes transporter level, the second and third terms 
on the right side of Eq. (UJ) can be rewritten as, 
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TABLE I: Characteristics of metabolic reactions and their £a, £l, and £jv when S = 10 . Marked * if a water molecule acts 
as a cosubstrate. Data for y/ {(AE Z ) 2 ) Z / (E z ) z from Ref. @], for Km from Ref. H, and for (S z ) z from Ref. []]|. 



Enzyme 


Substrate Product 


sj ((AE Z ) 2 ) Z /(E Z ) Z Km (mM) (S z ) z (mM) 


U 


& 




YEL042W* 


GDP GMP, phosphate 


0.235 0.1 0.39 


-0.271 


-0.260 


-0.278 


YJL005W 


ATP cAMP, diphosphate 


0.267 1.6 2.52 


-0.183 


-0.237 


-0.322 


YPL111W* 


arginine ornithine, urea 


0.241 15.7 50 


-0.244 


-0.261 


-0.291 



Therefore, the right side of Eq. ([7]) includes only the 
molecular levels measurable by currently available exper- 
imental techniques, and can be used to estimate enzyme 
kinetic parameters for single cells. A more complete de- 
scription of this idea can be found in the accompanying 
supplemental material 



11] 



Even in cases where the primary interest is the kinet- 
ics at a population level rather than at a single-cell level, 
does it matter to consider the effect of cellular variabil- 
ity analyzed so far? Since a specific value of £ can be 
changed responding to different conditions rather than 
kept constant, {Ki} will be changed also. This fact indi- 
cates that even a functional form of reaction kinetics at a 
population level will be distorted from its basal form. For 
example, the experiments with S. cerevisiae [3] suggest 
the relationship ((AE Z ) 2 ) Z /(E Z ) 2 Z « C 1 /{E Z ) Z +C 2 where 
C\ and Ci are constants (C\ « 1200 molecules/cell). By 
assuming that ({Av z ) 2 ) z is sufficiently small, Eqs. {I]) 
and lead to 
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(11) 

This result clearly violates the original form of the MM 
equation, since the MM equation is modified with the 
introduction of the second and third terms in (• • • ) of Eq. 
(fTTj). Although "effective" values of kinetic parameters 
fitted to the original MM equation might work for narrow 
ranges of (E z ) z and (S z ) z that are used for the fitting, 
it will still be unavoidable to witness the breakdown of 
the MM kinetics like Eq. (ITT1) in the face of significant 
changes of (E z ) z and (S z ) z [12j conveyed by severe intra- 
or extracellular condition changes. 

So far, we have used the MM equation for 
f(E z ,{S JZ }, {Ki}) assuming negligible molecular fluctu- 
ations under given E z and v z . A recent study has sug- 
gested that an active transport mechanism of substrates 
as well as the presence of competitive enzyme inhibitors 
may manifest the effect of such fluctuations [l3| , and this 
effect can be incorporated in our study through the use of 
the corresponding f(E z , {S 1Z }, {Ki}) instead of the MM 
equation. 

In this Letter, we demonstrate that reaction kinetics 
of a cell population can be systematically deviated from 
that of single cells by inevitable and significant variations 
in enzyme abundance. This result would not be only re- 



stricted to the case of MM kinetics focused on here; more 
sophisticated kinetic equations than the MM equation 
would also face such deviations. Our findings indicate 
that widely spread discrepancies between in vivo and in 
vitro kinetics might be attributed at least in part to cellu- 
lar variability, because previously known in vivo kinetic 
parameters have been mostly obtained from population- 
level experiments. We expect that the ultimate devel- 
opment of single-cell metabolomic analysis will greatly 
facilitate a precise determination of biochemical kinet- 
ics. In particular, such single-cell-level analysis can be 
applied judiciously to key parts of biochemical pathways 
of which operation is highly sensitive to their kinetic pa- 
rameters. 
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